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We compare the properties of transmission across one-dimensional finite samples which are 
associated with two types of "quantum diffusion", one related to a classical chaotic dynamics, the 
other to a multifractal energy spectrum. We numerically investigate models exhibiting one or both of 
these features, and we find in all cases an inverse power law dependence of the average transmission 
on the sample length. Although in all the considered cases the quadratic spread of wave packets 
increases linearly (or very close to linearly) in time for both types of dynamics, a proper Ohmic 
dependence is always observed only in the case of quasi-classical diffusion. The analysis of the 
statistics of transmission fluctuations in the case of a fractal spectrum exposes some new features, 
which mark further differences from ordinary diffusion, and enforce the conclusion that the two 
types of transmission are intrinsically different. 
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I. INTRODUCTION 

The onset of chaos in classical single-particle dynamics 
is often accompanied by deterministic diffusion, which is 
a type of stochastic motion which mimics to some extent 
a random walk, even in the absence of external random 
agents. Typical examples are the motion of a particle 
elastically bouncing in an array of fixed scatterers, and 
the diffusive energy absorption occurring in nonlinear pe- 
riodically driven systems. In both cases, characteristic 
transport properties appear, that can be more or less ac- 
curately described by an equation of the Fokker-Planck 
type. Whether similar transport properties persist on the 
quantum level is an important question, which lies on the 
borderline between Quantum Chaology and Solid State 
physics, and has received much attention on both sides. 
A most remarkable result in this context is that chaotic 
diffusion is inhibited by quantization. A classical one- 
dimensional diffusion is quantally reproduced only on a 
finite time scale, and is in the long run stopped by an 
interference effect very similar to Anderson localization. 
Equivalently, if the transmission of particles across a fi- 
nite sample of a quasi-one-dimensional disordered solid 
is considered, then the inverse dependence on the sam- 
ple length characteristic of a classical diffusive transport 
(Ohm's law) is quantally observed only on length scales 
which are small in comparison with the Anderson local- 
ization length. 

It was realized that, in a quasi-one dimensional case, an 
unbounded diffusive (or anomalously diffusive) spreading 
of wave packets is often associated with a singular con- 
tinuous spectrum of the Hamiltonian (or the Floquet op- 
erator in the case of periodically driven systems). This 
sort of a spectrum, regarded as a mathematical curiosity 
a long time, appears quite frequently in incommensurate 
structures, such as involved e.g. by quasi-crystals and by 
electron dynamics in crystals in magnetic fields; in such 
cases fractal singular spectra of zero Lebesgue measure 
have been found. In the light of the above findings, and 
insofar as the quasi-one dimensional case is concerned, a 
quantum "diffusive" regime, i.e., one marked by a linear 
(or close to linear) increase in time of the mean square 
displacement, can exist either due to an underlying clas- 
sical chaotic diffusion, developing on small scales in which 
localization effects are not yet apparent, or due to a frac- 
tal structure of the spectrum; no scale limitations being 
necessary in the latter case. 

Prototype systems for these two situations are the 
Kicked Rotator (KRM) and the Harper (HM) model, re- 
spectively. Each of these two models exhibits one type of 
quantum diffusion, but not the other; therefore, the two 
phenomena are in principle independent of each other. 
Still, the question remains, whether the "fractal" un- 
bounded diffusion may offer a quantum counterpart for 
the classical unbounded propagation in cases when both 
are present. A particular aspect of this question is what 
kind of transmission properties are associated with the 



"fractal" diffusion, and how do they match with the 
ohmic transport typical of classical diffusion. 

We have studied this question using the two mentioned 
models plus a third one, which in a sense interpolates be- 
tween them. This is the Kicked Harper Model (KHM), 
which displays at once a chaotic classical limit marked by 
deterministic diffusion (which is missing in the HM), and 
a "critical" quantum regime with a multifractal quasi- 
energy spectrum (which is missing in the KRM). We have 
compared the properties of quantum transmission across 
finite samples which are observed in the three cases. The 
analogies and differences we have observed supply useful 
indications for answering the questions raised before. For 
example, according to our results, the Ohm law seems to 
be rather exceptional inside the class of quantum systems 
endowed with "fractal" quantum diffusion. The finger- 
print of this macroscopic law would be then represented 
by the concurrent presence of a quasi-classical diffusion. 

Other differences emerge on analyzing transmission 
fluctuations, and lead to the conclusion that the two 
types of diffusion are marked by essentially different sta- 
tistical properties, and are in no apparent relation to each 
other. 

The plan of the paper is the following: in Sec. [n] wc 
review the mathematical formalism which allows for a 
general treatment of the kicked dynamics as a scattering 
In Sec. 



process. 
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we present the old (A) and the new 
(B,C) results about the transport properties respectively 
in the KRM, in the KHM and in the HM. The conclusive 
sec. IV is devoted to the analysis of the transmission 
fluctuations. 



II. DISCRETE-TIME SCATTERING PROCESSES 

We are interested in one-dimensional scattering pro- 
cesses in which free particles impinge on a " sample" of a 
finite length, whereby they are either reflected or trans- 
mitted. The particle dynamics inside the sample will 
have a nontrivial diffusive character, both in the clas- 
sical and in the quantum case. Such a situation can- 
not of course occur in strictly one-dimensional, conserva- 
tive Hamiltonian cases, but it can be realized with rel- 
ative ease if one considers a discrete-time dynamics in- 
stead of a continuous-time one - an artifact which has 
proven extremely useful in Nonlinear Dynamics. We 
shall presently outline the general formalism needed to 
describe a discrete-time scattering process. 

We consider the discrete time dynamics of a particle on 
a line, with coordinate q and momentum p. Throughout 
this paper we shall use dimensionless variables. The clas- 
sical dynamics of the particle is described by the iteration 
of an area-preserving map, of the form: 
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where the generating functions S, V will have different 
forms, depending on the considered model. In all cases, 
however, the function S will be even and periodic in 
the momentum p, with period 2k. Outside the region 
Qo < Q < <Zo + L the motion is uniform at con stant 
speed —S'(p) and momentum p. Thus eqs. (2.1) de- 
fine a scattering problem, in which free particles coming 
from infinity impinge on the scattering region, wherein 
they undergo "collisions" described by the "potential" 
V(q), until they re-exit from either end. With appropri- 
ate choices of the generating functions S, V, the dynamics 
inside the sample will have a chaotic character, leading 
to fast randomization in p and to diffusion in q. 

In this case, the scattering trajectories sensitively de- 
pend on various parameters, including qo, which defines 
the position of the scatterer on the line; both transmis- 
sion and reflection are affected in an unpredictable way 
by small changes of qo, because qo, together with the in- 
coming momentum, defines the initial conditions for the 
chaotic map whose iteration defines the dynamics of the 
particle inside the scatterer. 

For this reason, we can generate statistical ensembles 
by placing the scatterer at different positions; different 
locations of the scatterer of length L will define different 
statistical samples in the ensemble. This is in clear anal- 
ogy with the problem of electronic transmission through 
disordered solids. In that case different samples corre- 
spond to different realizations of disor der . 

The Quantization of the dynamics (2.1) is obtained in 
a standard way by iterating the quantum one-step prop- 
agator: 



(2.2) 



where we have introduced the adimensional Planck con- 
stant h, and V(q) is taken ^ only for q inside the 
scattering region. 

Due to periodicity in p, the Bloch theorem applies, and 
the evolution ( ^2 ) preserves the " quasi-position" a = q 
The value a of the quasi-position is a constant of the 
motion, therefore its value can be arbitrarily fixed. The 
quantum dynamics at fixed a can be studied, by using 
the quantization rule: 

d 

q a = -iTi— + a (2.3) 
op 

defined on 2-7r-periodic wave functions. Quantized in this 
way, the position assumes discrete values nh + a, (n inte- 
ger), that is, the quantum motion at given quasi-position 
takes place on a discrete one-dimensional lattice. We 
shall denote |n) the eigenstate corresponding to the n— th 
site of this lattice. The quantum evolution at fixed a is 
given by a unitary oper ator U a obtained by writing q a 
instead of q in eqn.( |2.2| ). Whereas U acts in the space 
of square-summable functions on the line, U a acts in the 
space of square-summable functions on the lattice. 

The unitary operator U a has eigenvalues e %x l h , and 
A (mod 2kK) is known as the quasi-energy. The quasi- 
energy is obviously conserved under the evolution. 



The quantum scattering theory stems, as usual, from a 
comparison of the long- time asymptotics of the dynamics 
generated by ( |2.2| ) with the free dynamics generated by 
Uq = exp(iS(p)/h). On mathematical grounds, U a is a 
finite-rank perturbation of Uq. This entails the existence 
of the M0ller wave operators ti T = lim^-too lPUq and 
thus ensures the well-posedness of the scattering prob- 
lem. In order to find the Scattering Matrix it is conve- 
nient to resort to time-independent scattering theory, and 
to seek the scattering eigenfunctions of U corresponding 
to a given value A of the quasi-energy. The eigenfunctions 
of the free dynamics corresponding to the same quasi- en- 
ergy are: 
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where pj, (j = 1, ...N\) are the roots of the equation 
S(pj) = A (mod 2kK). 

In this way for each quasi-energy A there are N\ degen- 
erate eigenfunctions; every one of them corresponds to an 
open scattering channel at the given quasi-energy. The 
number N\ of such channels depends on the quasi-energy 
A and on the specific form of the function S. 

The scattering dynamics can be visualized as follows : 
plane waves, eigenfunctions of Uq, are scattered by the 
sample, thus giving rise to reflected and transmitted 
waves. The crucial feature of this process is that quasi- 
energy is conserved, that is, S(p) is changed by an inte- 
ger multiple of 2kH. The construction of the 5-matrix 
follows standard rules of the scattering theory, via the 
Lippmann-Schwinger equation : 
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relating free waves Uq with true, distorted waves u + . In 
the above formula G+ is the "free" Green function : 
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By inverting the system (2.5) one obtains the distorted 
waves u+ and the 5-matrix, whose elements are given 
by: 
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where j,k = l,...N\ and \uj\ — l/\dS(pj)/dp\ is the 
density of states. Given the 5-matrix, the transmission 
coefficient is obtained by summing the squared moduli of 
the 5-matrix elements related to transmission in a given 
direction of propagation. 

The S— matrix is sample-dependent because the ex- 
plicit form of U a depends on the choice of the sam- 
ple, that is, on qo, through the values of the function 
A(q) = exp(iV(q)/K) at q = rih + a, (qa < q < qa + L). 
The string of these ~ L/% values has a more or less dis- 
ordered character, depending on the choice of the func- 
tion V(q); in any case, for any nontrivial choice of the 
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latter, different samples correspond to different strings, 
and different strings correspond to different realizations 
of " disorder" . There is a connection between the depen- 
dence of the 5-matrix on the sample and its dependence 
on the quasi-position a, because changing the latter by 
h is equivalent to shifting the sample by one site along 
the q— lattice. This suggests that the above described 
process of averaging over disorder may be equivalent to 
averaging over all possible values of the quasi-position. 
This is certainly true when 2tt/Tl is an irrational num- 
ber and A(q) is a 27r-periodic function - a case that oc- 
cupies a central position in this paper. In that case, 
any quantity related to scattering is also a 2ir— periodic 
function of a, and its average over all samples at given 
a = aa coincides with the average over all the values of 
a = ao + nh 77jog?(27t), which are obtained from ao by 
iterating the shift a — > ah mod(2n). Indeed, since this 
shift is ergodic for 2tt/?i irrational, such "time-averages" 
are independent of the starting point a$ an d coincide 
with "phase-averages", i.e., averages over a . 



III. TRANSPORT PROPERTIES 



A. The Kicked Rotator 



With the choice S(p) — bcosp, V(q) — \q 2 the map 



( p.ip describing the dynamics inside the sample becomes 
the Standard Map, which has been widely studied both 
in its classical and in its quantum version; the latter is 
associated with a well-known dynamical model, known 
as the KRM. 

The quantum KRM is a paradigm of the quantum phe- 
nomena associated with classical chaotic diffusion, and a 
short review of its main properties will provide a proper 
frame for the following discussion. At br 3> 1 the clas- 
sical Standard Map is fully chaotic, and the q— motion 
is similar to a random walk, with a diffusion coefficient 
D — j3(b)b 2 /2 where (3(b) is a known bounded function of 
b . On length scales much bigger than b, the distribu- 
tion f(q, t) of a statistical ensemble of orbits is adequately 
described by the diffusive approximation, which assumes 
f(q,t) to satisfy the diffusion (Fokker-Planck) equation 

The properties of the quantum KRM (at fixed quasi- 
position) crucially depend on the arithmetic nature of 
h. If h is a rational multiple of 47r, then the evolution 
operator (2.2) has an absolutely continuous spectral com- 
ponent that produces a ballistic spread of wave-packets 
on the q— lattice. If K is a "good" irrational number, 
accumulated numerical evidence supports a pure-point 
spectrum, and exponential localization of wave-packets 
on the lattice, with a localization length £ ~ b 2 /h (the 
properties of the KRM are reviewed in B). A continuous 
spectral component is also known to exist when H/At: is 
a "bad" irrational, well approximated by rationals; how- 
ever, neither the exact type of approximation required, 



nor the nature of the corresponding continuous spectrum 
are known. Thus in all known cases the long-time nature 
of the quantum propagation is different from that of the 
classical motion; while the latter is diffusive, the former 
is either ballistic or localized. Since the Correspondence 
Principle entails that in the limit h — > the classical dif- 
fusive motion must be recovered, there must be a time 
scale to{K) separating the quasi-classical diffusive spread 
of wave packets from the purely quantum regime, (bal- 
listic or localized). It is known that tu ~ D/h, which is 
the time required to diffuse over the localization length 
I. 

The scattering problem associated with the KRM was 
investigated in detail in ||. In the classical case, for 
b <C L the diffusive approximation is valid, and the nu- 
merically computed transmission coefficient T is in ex- 
cellent agreement with the theoretical value: 
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(3.1) 



which is found by solving the diffusion equation in the 
sample of length L, subject to a ppropriate boundary con- 
ditions 0]. On comparing (3.1) with the kinetic formula 
T sa £o/L involving the mean free path (mfp) £q, we are 
led to £o <~ b. Notice that T is a sample-independent 
quantity. 

In the quantum detailed numerical analysis was 

carried out for the strongly incommensurate case. For 
any quasi-energy A a transmission coefficient T(A) was 
computed following the procedure outlined in Sec. ||. 
T(A) is a fluctuating, i.e. sample-dependent, quantity. 
The main facts about the statistics of this quantity that 
emerge from numerical computations are: 

(1) (T(A)) qo (the average of T(A) over q Q ) is but weakly 
dependent on A, and is therefore close to the total trans- 
mission (the average of T(A) over A), which is the quan- 
tum counterpart for the classical transmission (3.1). It is 
on f act a decreasing function of L which comes close to 
( |3.lD in a range of lengths £q <C L <C £, or, equivalently, 
dClj'C b 2 /h, which is called the "diffusive regime" (see 
Fig.l). 

At larger lengths, (T(X)) qo exponentially decreases 
with L, as exp(— 2L/£); this is the localized, or "insu- 
lator", regime. 

(2) In the diffusive regime, Var(T(X)) (the variance of 
T(A)) is roughly independent of length Q. 

Both (1) and (2) reproduce well-known aspects of elec- 
tronic transmission through quasi one dimensional disor- 
dered solids; in particular, the scale-independence of the 
magnitude of transmission fluctuations observed in the 
diffusive regime is the celebrated phenomenon of Univer- 
sal Conductance Fluctuations (UCF). 



B. The Kicked Harper Model 

The KHM is obtained by choosing S(p) = 
dcosp, V(q) — ccosq. It was introduced in |M|J; the 
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study of its quantum diffusive properties was started by 
Lima and Shepelyansky ||, and was subsequently con- 
tinued by a number of Authors [p^|-p^[ (for a review on 
the subject see [jL5[ ). It is classically chaotic and diffusive 
at d, c S> 1, and the statistical analysis outlined for the 
KRM is still valid (with appropriate modifications). 

The spectral properties of the quantum KHM at fixed 
quasi-position, like those of the quantum KRM, depend 
on the arithmetic nature of h. If the latter is commensu- 
rate with 2-7T, then the spectrum is absolutely continuous, 
and the propagation of wavepackets is ballistic; a contin- 
uous spectral component is also present for "not too in- 
commensurate" values of Ti. However, even for strongly 
incommensurate values of h the spectral properties de- 
pend in a highly nontrivial way on the value of the pa- 
rameters c, d. Whereas exponential localization is still 
present at d > c, the motion is ballistic at c » d; a 
phase- portrait, illustrating the type of motion in various 
region of the (c, d) plane, can be found in [jl2| . These 
nontrivial properties of the KHM are connected with its 
double periodicity (in q and in p as well) JTT| , a feature 
which the KHM shares with the HM, which is the object 
of the next section. 

In particular, at c = d the motion is always delocalized, 
and wave-packets spread according to a law ((Aq) 2 ) ~ t~* . 
The exponent 7 is close to 1, but it appears to depend on 
the value of c = d. This behaviour signals t he pre sence of 
a singular-continuous spectral component flqjl 7| and in 
fact a multifractal spectrum has been numerically found, 
with a Hausdorff dimension close to 0.5. 

We have investigated the question, how is this (pos- 
sibly anomalous) diffusion at c = d reflected in the 
transmission properties of the scattering model. In 
our numerical computations we have used the function 
2arctan(ecos(-))/e with e << I as an approximation for 
the function cos(-) appearing in the generating functions 
S, V. The dependence of the classical transmission on 
the sample length is shown by the dashed line in Fig. 2 
and by the continuous one in Fig. 3; as in the KRM, it 
follows the law (3.1) with the appropriate parameters. 

The behaviour of the quantum (T(X)) qo depends on 
whether A is chosen in the spectrum, or not. In strongly 
incommensurate cases, the quantum spectrum is singu- 
lar, i.e., has zero Lebesgue measure, and (T(X)) qo is 
exponentially small at large L for almost all values of 
the quasi-energy A. Therefore, the total transmission, 
too, is exponentially small at large L, at varia nce with 
the classical behaviour, which follows the law (3J) in- 
stead. Thus, as in the KRM, the quasi-classical ohmic 
behaviour appears only in a limited range of lengths, 
which is larger, the smaller h. In Fig. 3 the dependence 
of (lnT(X)) qa on length is shown, for a strongly incom- 
mensurate case Ti/2tt = 1/(6 + p), where p is the golden 
ratio, p = (v5 + l)/2, and for a value of A chosen at ran- 
dom. After a relatively small interval in which a rough 
agreement with the classical behaviour is observed, the 
transmission decays exponentially. On comparing Fig. 3 



with Fig.l we see that the ohmic region is here narrower 
than for the KRM, though the ratio D/M , which in 
the KRM defines the width of this region, is roughly the 
same; moreover, fluctuations are definitely bigger. 

In order to detect the influence of the quantum dif- 
fusion associated with the singular spectrum we have 
to take A in the spectrum. A convenient strategy for 
locating the spectrum, and one commonly used in the 
study of incommensurate systems, consists in approxi- 
mating the actual, incommensurate system with a pe- 
riodic one. From a continuous fraction expansion one 
obtains a sequence of rational approximants p n /q n to 
H/2-k = l/(6+p). The KHM with H/2-k = p n /q n is then a 
periodic approximation to the strongly incommensurate 
KHM with h/2-K = 1/(6 + p). The spectrum of the in- 
commensurate KHM is approximated by the spectrum of 
the periodic KHM; the latter consists of bands separated 
by gaps, and the bands are generated as follows. Due to 
periodicity in q, the quasi-momentum j3 is also conserved, 
besides the quasi-position a. At fixed 0, a the dynamics 
is described by a unitary matrix of rank q n p[ ; as is 
varied at fixed a, the eigenphases of this matrix sweep 
the bands in the spectrum of the evolution operator U a , 
for the given a and for the given periodic approximation. 

Due to absolute continuity of the spectrum, the long- 
time dynamics of the periodic KHM is ballistic, (Ag) 2 ~ 
i 2 . Nevertheless, over fixed, finite scales of length and 
time the incommensurate dynamics will be reproduced 
by the periodic ones, the better, the higher the chosen 
approximant. This is illustrated in Fig. 4, where the de- 
pendence on time of (Aq) 2 (t) is shown for the infinite 
KHM for two rational approximants; two lines showing 
the ~ t and the ~ t 2 behaviour are drawn for compari- 
son. (Aq) 2 (t) was computed as the expectation value of 
q 2 over the wave function obtained at time t by iterating 
the KHM propagator 

JJ — p (ic/h) cos(q a ) e (id/h) cos(p) 

with a = on the initial state |0). 

Therefore, in order to analyze the scattering in a fi- 
nite sample, we have chosen a suitably high approximant, 
and have located the spectrum by a combination of two 
methods: first, by computing the bands by direct diag- 
onalization of matrices, as mentioned above; second, by 
analyzing the behaviour of the transmission T(A) as a 
function of A. An example of the results obtained with 
this procedure is shown in Fig. 5, where a scan of T(A) 
over a range of A encompassing one of the bands obtained 
by diagonalization is shown. 

Statistical ensembles were generated by varying a in 
the interval (0, 2ir). 

We have then studied the dependence of (— lnT(A)) a 
on length L for A in the center of a band; the result 
is shown in Fig. 2. The long-time ballistic propagation 
(entailing a length-independent transmission) is not yet 
manifest in the inspected range of lengths, as confirmed 
by the simulation of the infinite KHM. Open circles are 
well fitted by the full line of equation : 
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-<MT)) C 



-1.47 + 1.14 • \n(hL) 



giving once more a conductance inversely proportional 
to a power of length with an exponent close to 1, that 
could not be more precisely determined from our data. 
Thus the diffusion (with an exponent also close to 1) 
associated with the multifractal spectrum of the KHM 
results in an Ohm-like dependence of transmission on L. 

However, this quantum diffusive behaviour of trans- 
mission is quite different from the classical diffusive be- 
haviour, represented by the dashed line, and must there- 
fore be considered as a purely quantum effect. A quasi- 
classical behaviour is to be expected only in a region of 
very small L, (hardly detectable with the parameter val- 
ues used for Fig. 2), like in Fig.l. 

An obvious question is whether the difference be- 
tween the classical and the quantum transmission curves 
couldn't be simply imputed to a difference between the 
classical and the quantum diffusion coefficients. There- 
fore we have computed the quantum diffusion coefficient 
for the case of Fig. 2 from a direct simulation of the infi- 
nite KHM; the result was in fact different from the clas- 
sical coefficient ( D Q ~ D c i/3), but still, on substituting 
this value in eqn. ( |3.l| ) , we could not reproduce the quan- 
tum curve. We have therefore to conclude that the quan- 
tum transport associated with the multifractal spectrum 
is intimately different from classical diffusive transport. 

Finally, we emphasize that on choosing different quasi- 
energies in the spectrum different transmission curves are 
obtained (compare for instance open and full circles in 
Fig. 2). Therefore, averaging over quasi-energies (for 
example over the bands) destroys the Ohmic dependence. 



C. The Harper model 

In order to put the above results into a broader con- 
text, we now describe some results about transmission 
in the HM. This model does not belong in the class of 
kicked models discussed in Sec. || and does not possess 
a chaotic classical limit: still, under appropriate condi- 
tions it exhibits a multifractal spectrum and quantum 
diffusion. 

The Harper equation, also called almost Mathieu equa- 
tion, is a 1-d tight-binding equation: 

(H h ip)n = 4>n+i + ipn-i + V cos(Qna + f3)lp n = Elj) n 

(3.2) 

where ?p n is the wave function at site n, a is the lat- 
tice size, (3 is a phase, E the energy and Vb the poten- 
tial strength. Its diffusive properties have been studied 
by a number of authors (for exhaustive reviews see Ref. 
fl8f and Ref. jl9]]), and we shall here mention only a 
few characteristics connected with its transport proper- 
ties. Aubry and Andre [^(J have proved that this model 
presents a metal-insulator transition for Vb — 2, when 
Qa/2ir is a "strong" irrational number (duality theory). 



This means that if Vb < 2 the spectrum is absolutely con- 
tinuous with extended eigenstates, while for V > 2 the 
spectrum is pure point with localized eigenstates. At the 
critical point Vb = 2 the spectrum is singular continuous 



with zero Lebesgue measure 22 



If Qa/2ir = p/q, p, q integers, the Harper model is peri- 
odic |18| ; if Qa/2n is a "bad" irrational (a Liouville num- 
ber) Avron and Simon pl[| have proved that the spectrum 
is singular continuous (instead of pure point) for Vb > 2. 

The dynamical implications of these findings were 
studied in (l^j2^1 by solving the time-dependent 
Schroedinger equation: 

It was found that at the critical point Vb = 2 the spread- 
ing of wave packets obeys Var(n) ~ t 1 with 7 very close 
to 1 which means diffusive motion. 

The behaviour of transmission at the critical point has 
been investigated by many authors. Sokoloff [M pointed 
out that this model has, at the critical point, for Qa/2n 
a "strong" irrational number and E = (which belongs 
to the spectrum), very unusual transmission properties, 
namely wide fluctuations of the transmission coefficient 
as a function of the length (up to three order of magni- 
tude for variations of the length of order 1%). 

At the same time it was remarked that the 

transmission is an highly oscillatory function of the phase 
(3 too. 

An average resistance 



Rav(L) 



i=i 



where R(i) is the resistance of the sample of length i, 
was then introduced in Ref. Q and used by many au- 
thors, essentially as a means for distinguishing localized 
states from extended ones |p^ , p7| -|3lf . The average resis- 
tance was found to oscillate with the sample length for 
extended states and to increase exponentially for local- 
ized ones but, to the best of our knowledge, there were 
no claims about its behaviour as a function of length at 
the critical point. In Ref. |5(j) a Renormalization Group 
analysis was used to find a scaling relation for the average 
resistance near the critical point; the average resistance 
was found to diverge with the length as 2 — Vb — > + . In a 
closely related system, the Fibonacci Model, Sutherland 
and Kohmoto P9| have found a power law growth for the 
resistance. 

We have found that, if the average over the phase (3 
is taken, an Ohmic dependence is observed. To calculate 
the transmission coefficient we used the well known trans- 
fer matrix method (see for instance Ref. pi}|). As in the 
previously considered models, the scattering problem is 
one for waves propagating on an infinite 1-d lattice. Free 
waves have the form exp(ifcna) and the corresponding 
dispersion law is E = 2cos(fca). The potential Vb is dif- 
ferent from zero only for n — 1 , . . . , L — 1 . The reflection 
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and transmission coefficients are then found connecting 
the right-hand wave-function tp n = e ~ lkna -f re lkna f or 
n > L with the left-hand wave-function ip n = te~ lkna for 
n < 1. The transmission coefficient is given by T = \t\ 2 . 

We have fixed E — 0, (which belongs to the spectrum 
of the infinite HM because of symmetry), Qa/2ir = p = 
(\/E - l)/2, and F = 2 (critical point). In Fig.6 we show 
the average (— ln(T)) versus ln(L). The best linear fit 
is the dashed line, with slope 1.02 ± 0.03 and intercept 
0.5 ± 0.4. Thus we have an average "ohmic" law, with 
oscillations superimposed, which are associated with the 
rational approximants to the golden mean. In the same 
Fig.6 we also plot data related to three different ratio- 
nal approximants (the 7 th , 9 th and 13 th ) to the golden 
mean. These data follow the continuous line up to a 
length L ~ q (where p/q is the rational approximant); at 
larger lengths they saturate, i.e., transmission becomes 
independent of length, as expected of a pure periodic 
system. 

As a function of (3 , T shows erratic fluctuations, and 
the distribution of m(T) is well approximated by a log- 
normal distribution (for L > 1000), see Fig. 7. The log- 
normal distribution of the transmission coefficient is a 
nontrivial occurrence, because the spectrum is singular 
continuous and not pure point (the localized phase is 
known to exhibit a lognormal transmission distribution) . 

Averaging over (3 the logarithm of T rather than T 
itself was quite effective in significantly reducing fluctu- 
ations. Huge fluctuations of (T) p are indeed observed 
when the sample length is equal to a multiple of the de- 
nominator of the rational approximants to Qa/2ir, see for 
instance the full line on Fig. 8. This effect is particularly 
evident at small length, where only a few "resonances" 
are present. 

On checking the stability of these results against 
changes of the energy E and of the modulation Qa, two 
relevant aspects emerge. On one hand, no significant 
changes appear if the energy is changed inside the spec- 
trum, but far from the band edges, (this was done as in 
section 111(B)). On the other, a different irrational mod- 
ulation (quadratic or transcendental numbers) yields a 
different power law, as can be inferred from the Table 
H For brevity we call such a dependence a "generalized 
Ohm's law", and the corresponding exponent an "ohmic 
exponent (OE)". 

The OE in the third column of Table I is obtained from 
the best linear fit of (— ln(T))p vs ln(L); the error in the 
last digit is indicated in brackets, and the continuous 
fraction expansion is given in the 2nd column. 

The degree of irrationality of Qa/2Tt, rather than its 
absolute value, appears to determine the exponent of 
the algebraic decay of the transmission. In fact two 
different numbers only differing in the first terms of 
their continued fraction expansion (the golden mean and 
[100, 1, 1,1,.. .]) produce practically the same exponent. 
The fact that in Table I the exponent 1 (or very close 
to 1) occurs only for these two numbers raises the ques- 



tion, whether such an exponent is connected to a noble 
irrational modulation. Different choices of Qa, even in 
the class of the irrational algebraic numbers, produces a 
different OE. 

The last point is crucial, because, as reported in [p2[ , 
the dynamical exponent ruling the spreading in time of 
wave packets in the infinite Harper model is practically 
independent of Qa (provided it is sufficiently irrational), 
and so appear to be the scaling properties of the spectrum 

s 

This means that the OE is not directly connected with 
the spectral properties, and must therefore be determined 
by the structure of the cigenfunctions. Scaling proper- 
ties for the latter have been found |3(|, which however 
do not hold for generic irrational (e.g., transcendental) 
modulation, and cannot therefore by themselves explain 
the generalized Ohm's law, which is observed even for 
such numbers (bottom of Table I) . 



IV. TRANSMISSION FLUCTUATIONS 

In the previous sections we have studied the implica- 
tions that a diffusive spread of wave-packets, whether of a 
quasi-classical origin or due to a fractal spectrum, has on 
the average transmission. In this section we shall instead 
analyze the fluctuations of the transmission. 

In the theory of transmission through disordered solids, 
the quantum diffusive regime is characterized by Univer- 
sal Conductance Fluctuations, i.e., by transmission fluc- 
tuations of scale-independent magnitude |34| B7j PS ^TJ . 
On defining the dimensionless conductance as G = \T , 
where x is the number of scattering channels, it was found 
that 
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for a time-reversal invariant transmission process 
p7| , ^4 36 5l]] through a quasi-one dimensional sample of 
length L, with £q -C L -C I. Such universal fluctuations 
are considered a distinctive mark of a proper quantum 
diffusive regime, besides the validity of Ohm's law (apart 
from weak-localization corrections) . 

As mentioned in Sec. 0(A), the KRM exhibits UCF 
||||. In Figs. 9, 10 we present results for the other two 
models examined in this paper which do not share this 
property. Data for the KHM (Fig. 9) are rather ambigu- 
ous. In fact, even if the points are distributed along the 
line 2/15 (their average value is 0.12) they appear at the 
same time slowly decreasing with length. This behaviour 
is especially clear with the HM (Fig. 10), , where a power 
law dependence on length is found: Var(T) ~ L~ 5 with 
S = 0.38 ± 0.04 (in this case conductance and transmis- 
sion coefficient arc the same since only the "in" and "out" 
scattering channels are active). 

One more feature preventing fluctuations from being 
length-independent is the presence of sharp variations of 
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Var(T) as a function of the sample length L, when L 
is some multiple of a denominator occurring in rational 
approximants of Qa/2n. This behaviour does not ap- 
pear to be due to insufficient statistics, and is shown by 
a dashed line in Fig. 8. Therefore, although the range of 
lengths inspected in Figs. 9-10 is "ohmic" for transmis- 
sion, transmission fluctuations are not universal, because 
Var(T) is length-dependent. UCFs appear then to be 
typical of quasi-classical diffusion. 

Finally we remark that in the Harper model the trans- 
mission remains a non self-averaging quantity |5l] | , since 
Var(T) J (T) 2 diverges algebraically with L when L — ► oo. 



spectrum comes into play. 
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V. CONCLUSIONS 

The Ohmic dependence of transmission on the sample 
length is a typical feature of classical one-dimensional 
diffusive transport, which can be observed in strictly de- 
terministic classical systems in the presence of a chaotic 
dynamics. On the other hand, one-dimensional chaotic 
diffusion is typically suppressed, in the long run, by quan- 
tization, so that an ohmic behaviour can be quantally 
observed only on length scales which are small in com- 
parison to the quantum localization length. In some 
cases, however, localization effects leave room for an un- 
bounded pseudo-diffusive spread of wave packets, con- 
nected with a fractal structure of the energy (or quasi- 
energy) spectrum. On analyzing the transmission prop- 
erties associated with such a "fractal" pseudo-diffusion, 
which is observed in a model endowed with a chaotic 
classical limit, we have found evidence that this kind of 
quantum diffusion is a purely quantum effect, as unre- 
lated to quasi-classical chaotic diffusion as is quantum 
localization itself. On comparing the behaviour of the 
Kicked Harper model, which exhibits both types of dif- 
fusion, with the Kicked Rotator on one hand, and with 
the Harper model on the other, it becomes apparent that 
the proper ohmic regime is still restricted to a range of 
"small" sample lengths, in which the singular continuous 
nature of the spectrum is still unresolved (as it was the 
singular point-like spectrum of the Kicked Rotator), and 
the quantum motion is essentially a quasi-classical diffu- 
sion. Beyond that, the purely quantum pseudo-diffusion 
alone is at work, leading to a substantially different type 
of transmission, which is still power-law dependent on 
the length, but with an exponent which depends on the 
arithmetic nature of the Planck's constant, according to 
as yet unknown rules, and cannot therefore be meaning- 
fully carried to the classical limit. In addition, the statis- 
tics of transmission fluctuations are quite different from 
the quasi-classical case; in particular, their magnitude is 
not scale- independent, as it should in a proper metallic 
regime. In the language of the theory of quantum chaotic 
scattering, the crucial aspects of the Ericson regime of 
strongly overlapped scattering resonances appear to be 
lost as soon as the continuous (singular) nature of the 
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FIG. 1. Transmission versus the rescaled sample length (in 
units of the mean free path) for the KRM. Full symbols are 
quantum data for br = 10, b = 58, % = 1. Each point is 
obtained by averaging over an ensemble of 10 2 , 10 3 different 
samples. The quasi-energy A is fixed. Full line is the classical 
theoretical prediction (3.1). 



FIG. 2. Transmission versus the sample length for the 
KHM and c = d = 10, h/2ir = 144/1097. Full and open 
circles are relative to two quasi-energies in the middle of two 
different bands. Full line is the best fit for open circles. Each 
point is obtained by averaging over an ensemble of 10 2 differ- 
ent samples. Dashed line is the corresponding classical line 
from a kinetic equation like (3.1) with the classical diffusion 
coefficient. 

FIG. 3. Transmission versus the rescaled sample length (in 
units of the mean free path) for the KHM. Full symbols are 
quantum data for c = d = 10, h — 2n/(6 + p), p the golden 
mean. Each point is obtained by averaging over an ensemble 
of 500 different samples. The quasi-energy A is fixed. Full 
line is the classical theoretical prediction. 

FIG. 4. Energy growth in time for the KHM, for c = d — 10 
and rational H/2tt approximant to 1/(6 + p), p being the 
golden mean. Full line is for h/2n = 55/419, dots are for 
h/2-K = 144/1097. The quasi-position is fixed a = 0. For the 
sake of comparison t and t 2 lines are shown in dashed style. 

FIG. 5. Transmission coefficient as a function of 
quasi-energy for the KHM. Here c = d — 10, h/2ir = 5/38, 
L = 120. 



FIG. 6. Transmission versus length for the HM. Here 
E = 0, Vo = 2, Qa/2-K = p, p the golden mean. Each point is 
obtained by averaging over an ensemble of 10 3 random phases 
(3. Circles, squares and triangles are respectively for the 7 th , 
9 th , and the 13 th approximant to the golden mean. Dashed 
line is the best fit line with slope 1.02(3). 



FIG. 7. Distribution of the transmission coefficient for the 
HM, at fixed length L = 1000, by varying randomly [3 in [0, n], 
for E = 0, Vo = 2, Qa/2-K — p. The full line is the gaussian 
obtained from the standard best fitting procedure. 

FIG. 8. /3-averaged transmission coefficient (full line) and 
its variance (dashed line) as a function of the length L for the 
same data as Fig. 7 As vertical lines are shown multiples of the 
rational approximants to p: 13, 21, 34, 55, .... Averaged values 
were computed by taking 5000 random (5 in the interval [0, n]. 

FIG. 9. Var(G) versus the length for the KHM. Open and 
full circles are obtained by fixing two different quasi-energies 
and varying the quasi-position over 100 samples; triangles are 
obtained by varying the quasi-energy over 100 samples and 
fixing the quasi-position a — 0. Horizontal dashed line is the 
theoretical value 2/15. 
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FIG. 10. hn(Var(T)) versus the logarithm of the length for 
the HM. Here E = 0, V = 2, Qa/2n = p. Each point is 
obtained by averaging over 10 3 different samples. Dashed line 
is the best fit line with slope S = -0.38 ± 0.04. 



TABLE I. Generalized Ohm's exponent for various irra- 



tional modulations. 


Qa/2-K 


Expansion 


OE 


(V5-l)/2 


[1,1,1,1,...] 


1.02(3) 


(197- a/5)/19402 


[100,1,1,1,...] 


1.00(5) 


y/2- 1 


[2,2,2,2,...] 


1.35(4) 


a/3-1 


[1,2,1,2,...] 


1.15(8) 


(V3-l)/2 


[2,1,2,1,...] 


1.15(3) 


(V37-4)/3 


[1,2,3,1,2,3,...] 


1.28(4) 


e-2 


[1,2,1,1,4,...] 


0.75(3) 


7T-3 


[7, 15,1,292,1,...] 


0.38(3) 
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